Isoscaling bearing information on the nuclear caloric curve 
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Abstract 

We show that the qualitative behavior of the nuclear caloric curve can be inferred from the 
energy dependence of the isoscaling parameters. Since there are strong indications that the latter 
are not distorted by the secondary decay of primary hot fragments, in contrast to other observables, 
this suggests that valuable information on the nuclear caloric curve may be obtained through the 
analysis presented in this work. 
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The determination of the nuclear caloric curve is key to the understanding of the multi 



fragment emission observed in heavy-ion collisions at intermediate energies 

HQ 

. Indeed, 

equilibrium statistical calculations predict a continuous rise of the breakup temperature as 
a function of the excitation energy of the disassembling system if the multifragment emis- 
sion takes place at fixed density P, S, Q|. In this case, the pressure at the breakup stage 
increases monotonously as a function of the excitation energy. On the other hand, a wide 



plateau should be observed 



if the onset of multifragmentation occurs after a rapid rise 



of the pressure, which is accompanied by a steady decrease of the breakup density as the 



excitation energy increases 15] . Similar conclusions have also been obtained with a soluble 
thermodynamic model by Das et al. |2J. 

The experimental observation of the caloric curve has been intensively debated in the last 
vears fflflflfl since different r_en t s nave led to distinctly different conclusions 
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[6|, [10J, Hil, [12J, HJ, [15J, Ha]. One problem is due to the great difficulties in measuring 
the nuclear temperature as this quantity can only be inferred from fragment information 
measured long after the breakup stage. It has been shown that side feeding from the de- 
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excitation of the primordial hot fragments 
drawn from the different methods currently employed [18, 
of the breakup temperature. 



may bias empirical conclusions 



20 



22] in the determination 



In this context, the isotopic scalin g re cently observed in nuclear reactions, in a broad 



range of bombarding energies 



H 24,0,0, 



27, 



, is expected to be rather insensitive 
to effects associated with the secondary decay of the primary fragments More 
precisely, the ratio R21 between the multiplicity Yi(N,Z) of a fragment whose proton and 
neutron numbers are, respectively, Z and N, measured in two reactions (1 and 2) with 
different isospin, follows the relation: 

Y 2 (N, Z) 



R 
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Cexp (aN + (3Z) 



(1) 



Y!(N,Z 

where a and (3 are the isoscaling parameters and C is a normalization constant. This scaling 



27, 



281 ] . it has been observed in deep inelastic 



property is very robust. As is shown in refs. (24 , 
reactions, evaporation processes, besides nuclear multifragmentation. 

One explanation for this scaling behavior may be found in the grand-canonical ensemble, 
in which the multiplicity Y(N, Z) is given by: 

/ipZ + ji n N 



Y(N,Z)=Uz(T,V f )exp 



T 



(2) 
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where T is the breakup temperature, \i n and \i p respectively stand for the neutron and 
proton chemical potentials and 



(az = g A z^-A 3/2 exp 

Arp 



Baz — Jaz{T) 



(3) 



T 

In the above expression, A is the mass number, Baz and qaz are, respectively, the fragment's 
binding energy and spin degeneracy factor. The excitation energy of the fragment is taken 
into account by the internal free energy Jaz{T). The free volume Vf is a parameter of the 



calculation and A^ = y2Txh 2 /m^T 1 where is the nucleon mass. 

If the breakup takes place in the two reactions at the same temperature and density, the 
ratio involving (az cancels out and one finds that: 

(2) _ (1) (2) _ (1) 

a = — and (3 = ^ , (4) 

where the superscripts label reactions (1) and (2). Although the observed yields are affected 
by the decay of the primary fragments, the form of Eq. (0) is still expected to hold, even 
though (az will be given by a more complex expression. Since (az might be very similar 
in the two reactions, as long as the breakup temperatures are close enough, the isoscaling 
parameters should be safely obtained from the final yields [3, 2^] and, therefore, should 
bear reliable information on the breakup stage. 

Owing to the exponential relationship between the system mass and the chemical po- 
tential, // 2 ) — fjP^ should be a slowly varying function of the temperature. This is indeed 
confirmed by our calculations (see below), in agreement with the findings of ref. [26]. There- 
fore, the robustness of the isoscaling parameters may be used to investigate the qualitative 
behavior of the caloric curve. 

To address this point, we apply the isoscaling analysis to the decay products of reactions in 
which a thermally equilibrated source with excitation energy E* breaks up statistically. We 
consider a proton rich source, whose mass and atomic numbers are A = 168 and Z = 75, 
which corresponds to reaction (1). In the second reaction, we use a neutron rich source, 
with A = 186 and Z = 75. The improved version |3| of the statistical multifragmentation 
model (ISMM) developed in refs. ^,0,0], and described below, is used in the following 
calculations to simulate the decay of the excited source. Either implementation of this model 
is useful in the present study since, under different assumptions for the breakup density, it 
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predicts qualitatively different caloric curves 



Thus its predictions can provide the input 



to investigate the sensitivity of the isoscaling parameters on the qualitative shape of the 
caloric curve. 

In the ISMM model, partitions strictly consistent with the constraints: 

Ao = ^2N AZ A, Z = Y,N A zZ, (5) 

AZ AZ 

and 

+ E* = \^ + Y, N AZ E AZ (T, V) , (6) 
5 K ° AZ 

are imposed. In the above equations, N AZ denotes the multiplicity, in each generated par- 
tition, of fragments whose mass and atomic numbers are A and Z, Eq S ' is the ground state 
energy of the source, e represents the elementary charge, and Rq is the radius of a sphere 
with a volume V, corresponding to the breakup volume. The energy E AZ (T,V) contains 
contributions from the fragment's binding energy, excitation energy, translational motion, 
besides the remaining Coulomb terms which, through the Wigner-Seitz approximation, pro- 
vide the corrections to account for the Coulomb repulsion between the fragment s |30 j. A 
Monte Carlo sample of the possible fragmentation modes is carried out following |2]j. The 



breakup temperature is determined, for each partition, by solving Eq. 

The main differences from the ISMM U and the original SMM 3(J, |3jj are in the 
use of internal free energies built from empirical data on discrete states wherever available 

n 

and of experimental binding energies all over the mass table [221]. Careful extrapolations 
are carried out to mass regions where the information is not available in either case ^ 13 ■ 
Both quantities influence directly the determination of T in each fragmentation mode /, 
whereas the free energies also play an important role in the evaluation of the entropy Sf, 
which enters in the calculation of any physical observable AZ : 



f AZ exp 


T,{AZ} f n azS az 


E/ ex P 


^{Az} f NazSaz 





(Oaz) = — r 1 (7) 



Due to the constraints imposed on each partition, a physical observable AZ fluctuates from 
one fragmentation mode to the other. 

The average breakup temperature may be obtained through the above expression and it 
is shown in Fig. das a function of the excitation energy for the two sources we consider. 
In order to suppress statistical fluctuations, a billion events have been generated for each 
excitation energy. A monotonous increase of the temperature is observed in the upper 




E / A (MeV) 

FIG. 1: Caloric curve predicted by the ISMM for the neutron (triangles) and proton (squares) rich 
sources. In the upper panel it is assumed that the breakup occurs at a fixed density whereas the 
breakup volume is allowed to change in the results displayed in the bottom panel. For details see 
text. 



panel, in which case the breakup volume of the system is kept constant and is 6 times larger 
than its value at normal density. One also observes that the average temperatures are very 
similar for both sources, which justifies the assumption that the temperature dependent 
terms in (az cancel out in Eq. As stated before, and reported in ref. the situation is 
qualitatively different as one allows the breakup volume to be multiplicity dependent and a 
plateau is observed approximately between 3.5 and 7 MeV/A. The breakup temperatures in 
both systems are quite close, except at the region around the onset of multifragmentation, 
E*/A ~ 3.0 MeV, where small differences may be noticed. 

In contrast to an earlier study also using the SMM Q, in which a larger isospin depen- 
dence in the plateau region was predicted, such effect is rather reduced in our calculations. 
To understand this aspect, we rewrite Eq. (0) in the following way: 



= l T (M-l) + Y J N A E\(T). 

Z A 

(8) 

where N A denotes the multiplicity of fragments with mass A in the fragmentation mode 
considered and the total fragment multiplicity is represented by M. The quantity R ce u 



E* - 



B 



AqZ 



N AZ B AZ 



AZ 



3Zy ^ 3ZV 



5 R 



AZ 
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cell 
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stands for the radius of the cell in which the fragment is embedded as one subdivides the 
system in order to apply the Wigner-Seitz approximation [soj ]. 

The first term on the right hand side of Eq. (jSJ) corresponds to the fragments' kinetic 
energy (in the center of mass reference frame of the total system) and E* A (T) represents the 
internal excitation energy of a fragment with mass A at temperature T. We suppressed any 
Z dependence in E* A (T) since it has been introduced in SMM only in ref. jjj. It entered in 
previous calculations only in the case of A = 4 nuclei because all light fragments (A < 5) were 
assumed to have no internal degrees of freedom, except for the alpha particles. Therefore, 
all the dependence on isospin is confined to the left hand side of this equation. 

An increase of Z /A enhances the Coulomb term between brackets in Eq. (JSJ) and, as a 
consequence, tends to lower the temperature, which appears only on the right hand side of 
the above expression. This is the effect observed in ref. j^. 

Nevertheless, the binding energies play a very important role in the balance of the equa- 
tion and, consequently, on the temperature extracted from it. It is shown in ref. 2^ that 
total binding energies calculated through simple liquid drop mass formulae, such as that 
used in ref. j23j, deviate appreciably from empirical values. The differences reported in ref. 
(2^] can be as large as 40 MeV for heavy nuclei and are, on the average, around 10 MeV for 
light fragments. Therefore, Ba z — J2az NazBaz calculated in [33[ has systematic errors 
which further increase the isospin dependence of the temperature. Our results, which are 
obtained with empirical binding energies and careful extrapolations to mass regions where 
this information is not available, as described in 0, exhibit a much weaker isospin 
dependence. 

We now investigate the extent to which the isoscaling parameters carry informaton on the 
caloric curve and show, in Fig. |2l 1/aasa function of the excitation energy. The parameters 
are obtained by fitting Eq. to the primary yields predicted by the microcanonical ISMM 
calculations, with the same isotopes, Z = 1, 2,..., 8, used in ref. 1251. We have checked 

nn n 

that, in agreement with [25J, |29j], the changes in 1/a due to the secondary decay P( of the 
primordial fragments are of the order of 5% at E* / A = 3, 6, and 9 MeV and, therefore, do 
not change our conclusions. As in the previous plot, the results shown in the upper panel 
correspond to a fixed breakup density whereas it is allowed to change in the lower part of the 
picture. The results reveal distinct qualitative behaviors in each case. More specifically, the 
reciprocal of a follows approximately a straight line for almost the full range of excitation 
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FIG. 2: Reciprocal of the isoscaling parameter a as a function of the excitation energy. Constant 
breakup volume is assumed in the upper panel, while in the lower one the results were obtained 
with variable breakup volume. The lines illustrate the slopes in each region. For details see text. 



energies considered if the breakup density is kept fixed. A clear change of slope, before 
and after the plateau region, appears if the breakup volume is multiplicity dependent. The 
plateau observed in the caloric curve is not so apparent in the 1/a plot because, although 
/i^ 2 ) — //W varies slowly compared to T, it somewhat distorts the curve. The relevant point 
here is that one should observe different qualitative behaviors in the 1/a plot according to the 
characteristics of the caloric curve. Similar conclusions are obtained with the f3 parameter, 
but we concentrate on a because, as shown in refs. 0, I29 1 . f3 might change more than 
a after secondary decay of the primary fragments, probably due to Coulomb effects. The 
distinction between the two scenarios allows one to ascertain the existence of the plateau 

n 

in the caloric curve. It is worth mentioning that the experimental results presented in 26J 
seem to favour the monotonous increase of the breakup temperature. However, since in this 
reference 1/a is plotted as a function of the bombarding energy, further analysis is needed 
to draw more precise conclusions. 

In order to investigate the dependence of the results on the statistical ensemble employed, 
we also calculate the primary yields, for the two sources considered here, in the framework of 
the grand-canonical approach. To prevent artificial deviations, we use the same ingredients of 
the micro-canonical case, such as binding energies, internal free energies, and spin degeneracy 
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FIG. 3: Same as fig. EJ The circles represent the predictions of the micro-canonical approach, 
whereas the results obtained with the grand-canonical ensemble are depicted by the squares. For 
details see text. 

factors. The free volume and the breakup temperature, which enter into eqs. (j2J) and (|3*j). 
are obtained from the micro-canonical calculation for each excitation energy, instead of 
considering them as free parameters as is usually assumed. Given these two quantities, the 
chemical potentials are calculated by imposing that: 



A = ^Y(N,Z)A and z o = ^Y(N,Z)Z , (9) 

NZ NZ 

where Y(N, Z) is computed through eqs. (J2J) and ©. The chemical potentials are found by 
minimizing the difference between the left and the right hand sides of the above equations. 
We impose that the constraints are fulfilled with a precision better than 3 digits. Once the 
chemical potentials are obtained, the primary yields Y(N, Z) may be evaluated through eqs. 
© and ©. 

The predictions of the two ensembles are compared in Fig. |3J where the results of the 
grand-canonical calculations are depicted by the open squares whereas the full circles repre- 
sent those obtained with the micro-canonical ensemble. The isoscaling parameters are ob- 
tained by fitting Eq. using the grand-canonical yields, as is done in the micro-canonical 
calculation. It may be noticed that, although the absolute values differ from one statistical 
approach to the other, both exhibit the same qualitative behavior. We have checked that a 



S 



and (3 obtained through Eq. (jlj) agree fairly well with the grand-canonical values displayed 
in this picture. Small deviations are observed only in the regions in which the difference 
between the breakup temperatures of the sources are non-negligible. Besides showing the 
consistency of the isoscaling assumptions, this confirms the expectation that small differ- 
ences between the average temperatures of the two sources should not appreciably affect the 
scaling properties. 

The difference between the predictions of the two statistical approaches may be attributed 
to the strong mass, charge, and energy constraints imposed in the micro-canonical calcula- 
ions in each partition. These restrictions lead to a fairly broad temperature distribution 
34J], corresponding to different fragmentation modes. In contrast, the temperature is kept 
constant in the grand-canonical ensemble whereas the charge, mass and energy are fixed 
only on the average and not event by event. However, the qualitative agreement between 
the two calculations is a very positive aspect since it is extremely difficult to select events 
whose decaying sources strictly obey these constraints in experiments. Therefore our anal- 
ysis seem to indicate that, at least, the qualitative behavior of the caloric curve may be 
studied experimentally. 

In conclusion, we suggest that the existence of the plateau in the nuclear caloric curve 
may be better investigated experimentally through the isoscaling analysis. Clear deviations 
from the linear behavior are expected to be found in the 1/a vs E*/A curve if the plateau 
exists. Conversely, if this plateau does not exist, the 1/a vs E*/A curve should follow a 
straight line over a wide excitation energy domain. Our results indicate that the deter- 
mination of the qualitative shape of the caloric curve can be done more reliably using the 
isoscaling analysis than through measurements of the temperature from the multiplicities of 
the detected fragments. 
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